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, ^ , . Abstract 

New vortex solutions to the Landau-Ginzburg equations are described. These configurations, 
which extend the well known Abrikosov and giant magnetic vortex ones, consist of a suc- 
cession of ring-like supercurrent vortices organised in a concentric pattern, possibly bound 
10 i to a giant magnetic vortex then lying at their center. The dynamical and thermodynamic 

\l ■ stability of these annular vortices is an important open issue on which hinges the direct 

I experimental observation of such configurations. Nevertheless, annular vortices should affect 

' indirectly specific dynamic properties of mesoscopic superconducting devices amenable to 

I physical observation. 
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1. Introduction. Ever since they were introduced, the Landau-Ginzburg (LG) equations 
have remained a powerful and insightful tool for the exploration and understanding of supercon- 
ducting and other collective quantum phenomena, in particular in low temperature supercon- 
ductors. Among a host of important results and developments which were achieved through the 
use of these equations, the existence of Abrikosov magnetic vortex solutions [||] still stands out as 
a most remarkable prediction transcending the boundaries of condensed matter physics by itself, 
whose physical significance and applications still have many riches in store even forty years on, 
especially with the recent advent of mesoscopic and nanoscopic superconducting devices. 

In this letter, we wish to report on new solutions to the LG equations which share some 
properties with usual Abrikosov configurations, but are nevertheless of quite a different nature 
and are possible because of the boundary conditions implied by mesoscopic samples of finite size. 
Rather than having the order parameter vanish along a line within the superconducting material 
as occurs for Abrikosov and giant magnetic vortices, these solutions are characterized by having 
the order parameter vanish on one or more surfaces each with the topology of a cylinder and 
enclosing one another in a concentric-like pattern with more or less constant spacing, while the 
maximal number of such surfaces is set by the size of the mesoscopic sample relative to the 
coherence length ^ of the material^. Moreover, the volume defined by any pair of successive 
surfaces is threaded by a super cur rent running in a closed loop, parallel to and vanishing on 
these surfaces, with all these successive sheaths of currents having the same circular orientation. 
Finally, these configurations may also carry a non zero fluxoid quantum number, in which case 
they may be viewed as being bound to an ordinary giant vortex carrying that quantum of 
fiux and lying at the center of the set of surfaces, the associated linear locus of the vanishing 
order parameter then threading along the central axis of the innermost surface. Because of this 
particular topology, we refer to these solutions as "annular vortices"^ of order n and fluxoid 
L, n referring to the number of cylindrical surfaces on which the order parameter vanishes and 
L being the usual fluxoid quantum number. These solutions thus generalise the well known 
giant magnetic vortex ones which correspond to n = and an arbitrary number L of Abrikosov 
vortices or anti- vortices lying on top of one another [^. 

Since these configurations solve the LG equations, they define local extrema of the free 
energy, but it has not yet been established whether these are local minima rather than maxima 
or saddle points whatever the values for n and L. Irrespective of this issue of their dynamical 
stability in configuration space^, it is clear that for a fixed value of L, the free energy of these 
solutions increases with increasing n. Hence in practical physics terms, such configurations have 
only a finite thermodynamic lifetime which also still needs to be assessed. 

Consequently, at this stage, it is not clear whether these annular vortices could be ob- 
served in mesoscopic devices, and how low a temperature would be required for that purpose. 
Nevertheless, from the quantum mechanical point view, annular vortices — even if dynamically 
unstable — provide for new contributions to the density of states of such systems, and should 
therefore affect specific physical properties of mesoscopic devices. One such instance could well 
be the behaviour in an external magnetic field of mesoscopic disks and annuli||5|, whose maze 
of (L, n = 0) energy curves in their phase diagram is then filled by many more such curves for 
n 7^ 0, leading to a plethora of degeneracy crossings as a function of the applied field and thereby 
a dynamics of these devices possibly much richer still than that so far considered[^, ^. 

In fact, having in mind eventual applications such as a new type of quantum detector for 
elementary particles and atoms, or as the basic qubit device for a quantum computer, these 

few weeks after these solutions had been found, we received Ref.j^ in which similar ones are discussed in 
the context of quantum optics, based on coupled non linear equations analogous to the LG ones. 
■^The idea of a vortex being associated to that of a current running in a closed loop. 

^Note that the stability of the ordinary giant vortex solutions (L, n = 0) seems to have been settled only 
recently 1^. 
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annular vortex solutions were found precisely when studying mesoscopic samples with a cylin- 
drical and annular geometry, and assuming an infinite extent along their symmetry axis. In that 
case, in a plane perpendicular to that axis, the order parameter vanishes at the center when 
L 7^ 0, while a series of almost regularly spaced n concentric circles of vanishing order parameter 
appear when n 7^ (in the case of an annulus, the same description applies with the central hole 
simply removed). This is the situation presented in this letter, but it is clear that for samples 
whose boundaries are of arbitrary shape, the same classes of solutions exist, obtained through 
continuous deformations of the perfectly symmetric ones constructed here, thereby leading to 
the general picture for these annular vortices as described above. 

2. The LG equations. For time independent configurations in the presence of an external 
magnetic field Bext-, in S.I. units the free energy of the system writes asB 
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The notation used should be self-explanatory for most variables. The symbol stands for the 
volume of the superconducting sample, = 2Trh/\q\ — with q = —2e (e being the positron 
charge) — is the usual unit of quantum of flux, while A and ^ are the penetration and coherence 
lengths, respectively. However, we need to emphasize the fact that in this expression the order 
parameter iIj{x) is normalised to its constant value in the bulk of the material in the absence 
of any magnetic field, and that the free energy is defined such as to vanish exactly at the 
normal-superconducting transition. 

The above expression of the free energy makes it clear that it is the penetration length 
A which weighs the relative contributions of the magnetic field energy — due to deviations of 
B{x) from the applied field Bext — and of the condensate energy — due to deviations of |?/'(x)| 
from zero — , while it is the coherence length ^ which weighs the contributions to the condensate 
energy due to spatial inhomogeneities in tp{x) — through the covariantized gradient — relative to 
those due to deviations from the bulk value \^lJ\ = 1 — through the LG potential energy. This 
interplay of scales is made manifest in (|l]). 

Extremizing the free energy leads to the LG equation 
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K = A/^ being the LG parameter, coupled the Maxwell equation 

V X B = flQ Xm , 
with the electromagnetic current density given by 
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In particular, this latter relation implies the second London equation for the Meissner effect 
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Let us now particularize these equations to axially symmetric configurations, the external 
homogeneous field -Bext being also aligned with the symmetry axis iz- We shall specifically 
consider a sample either of cylindrical geometry with radius (and = 0) or of annular 
geometry with inner and outer radii and rf,, respectively. Finally, in the present study 
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we assume the extension of the sample in the zbe^ directions to be infinite. The following 
parametrisations then apply, 

B = B{r)e, , A = Air)e^ , ^oXm = J(r)e<^ , ^(r,(/>) = /(r)e^^('^) , p(r) = |V'(r,(/>)|2 = f{r) , 

(6) 

where (r, 0) define the polar coordinates in the plane perpendicular to the symmetry axis. 
In contradistinction to other works, an important point needs to be emphasized concerning 
the parametrisation used for the order parameter in terms of the two functions /(r) and 9{(j))- 
The function /(r) is of course real, but of a sign which may be positive or negative as the 
case arises, while the phase function 9{(j)) is assumed to be continuous. Usually, the function 
/(r) is explicitely assumed to be the square-root of the (normalised) condensate density p(r), 
/(r) = ^/p{r), whose sign is thus always positive, with the important consequence that whenever 
the order parameter vanishes such that the gradiant of does not vanish as well, on the "other 
side of the zero" there is necessarily a discontinuous jump of ivr in the phase factor ip which is 
then used to parametrise the order parameter through ip = y^e**^. There is of course no "other 
side of the zero" in the case of a giant vortex, but there is obviously one when ip vanishes on a 
surface within the sample^. The general parametrisation of ip{x) advocated here allows for the 
specific possibility that f{x) may change sign when crossing a zero of '0(f) = f{x)e'^(^\ all m a 
continuous fashion for both functions f{x) and 9{x) and their gradients. This point is obviously 
essential to our results. 

To proceed further, let us choose to normalise quantities in terms of the penetration length 
and the associated unit of magnetic field <I>o/(2vrA^), by introducing the dimensionless variables. 



The coupled LG-Maxwell equations then reduce to|8| 
1 
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while the magnetic field B{r) is determined by the second London equation 



h{u) = ^^g{u) . (9) 
u du 

Furthermore, a careful consideration of the boundary conditions (b.c.) of the problem, of the 
small u behaviour of solutions to these equations, of the flux threading the sample and of the 
relation (^) between ^0'^m and ip, shows that the remaining quantities are determined to be[P 

e{(f)) = -L(j) + eQ , ua{u)=g{u)+L , L = 0, ±1, ±2, . . . , (10) 

where is arbitrary and L is the usual fluxoid quantum number, thus corresponding to the 
order parameter ip{r,(f>) = /(r)e~*'^'^e*^" characteristic of giant vortex configurations for L 0. 
In particular, the total flux threading a disk of radius Ua < u < Uf, within the sample is 

q>{u) = $0 [giu) + L] = ^>o ua{u) , Ua < u < Ub . (11) 

The advantage of expressing the problem in these terms rather than those usual in the literature, 
is that only gauge invariant quantities still appear in the coupled equations to be solved, while 
the gauge variant ones, 9{(j)) and ua{u), are explicitely known already. 

''From that point of view, annular vortices are very much reminiscent of domain wall configurations in the 
spontaneously broken Yang-Mills theories of particle physics and cosmology, and could thus possibly be of interest 
in these other fields of physics as well, not to mention of course other coherent quantum systems in condensed 
matter physics. 
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The relevant b.c. still need to be specified. In the disk case surrounded by the vacuum or 
an insulating material, they are|^, |8[ 



at n = : g(n)|_„ = -L ; = if L = or /(n)|_„ = if L / , 

U, UiU, LLLL ll — 7/., 



(12) 



while in the annulus case, we have! 



^tu = u, : =6ext ; ^1 =0 . 



(13) 



Finally, it may be worth noting also that both these sets of b.c. as well as the coupled 
equations (^) for f{u) and g{u) derive from the variation of the following quantity (valid both 
for the disk with Ua = and for the annulus with Ua 7^ 0) 



" _ 1 {2tt>^\'^ F = F 

~ 2 2^ V *0 / ^norm , 
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where £ denotes the free energy E per unit length in the direction of the sample. 

3. Analytical insight. The resolution of the equations (|8|) requires a numerical approach, 
both because of the non linear character of these equations as well as the b.c. at finite values 
of r. Nevertheless, some analytical considerations may be developed, giving already the insight 
necessary to the construction of annular vortices and the actual demonstration of their existence. 

The pair of coupled second order equations (^) requires a total of four b.c. However, 
in order to determine a solution uniquely, these four conditions must be specified at the same 
boundary, thereby enabling the construction of the then unique solution by propagating these 
boundary values throughout the sample volume using the equations (^). In the present case, 
we do have a total of four such conditions, but split evenly between the two boundaries of the 
sample at Ua {ua = in the disk case) and U},- In order to construct solutions, one needs to 
add two free b.c. at either boundary, propagate the solution within the sample and then adjust 
the two free conditions in order to meet the two conditions specified at the other boundary. 
Consequently, this procedure may generate more than a single solution. 

To be specific, consider first the disk case and introduce the parametrisation 

f{u)=u\^\f{u) ; /o = /(n)u.o > = = ^(^)|"-o ■ (15) 

From the equations dH), it follows that f{u) and g{u) are both even functions of u, while the two 
boundary conditions at u = in (|l^) are then also satisfied. Thus, the parameters /o and qq 
correspond to the two additional free boundary values at « = in order to construct a unique 
configuration through (|^), with in particular related to the value 6ext of the external field. 
From a straightforward linearisation of the LG equation, one then finds exactly 



1 /2\l-^l 

^/(n) = \L\\ - , .y-; - , 



lim ^ — /(n) = |L|! ( - ) J\L\ii^u) , ^lim ^g{u) = -L , _^lim^&ext = , (16) 



where J^n{x) is the Bessel function of the first kind. In other words, precisely at the normal- 
superconducting transition in a vanishing external field, there is in fact lurking in the still waters 
of f{u) = an oscillatory pattern just waiting to emerge as soon as /o / 0, namely as soon as 
the superconducting phase is reached. Thus for /o arbitrary close to zero but not vanishing, by 
continuity there exists a solution to (^) and the b.c. ([l^ ) at u = for which f{u) passes through 
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a succession of zero values while alternating in sign. In the limit (/q = 0, go = 0), these values 
are the zeros of the Bessel function J\i^\{k,u). To a good approximation, their spacing in the 
sample is such that Ar ~ tt^, while the number of those lying within the radius of the sample 
defines an upper bound on the maximal number n^.^^{L) of annular vortices possible, differing 
from it by at most one unit since the b.c. ( |l2|) at n = Uf, still need to be solved, and whose value 
nmaxiL) thus decreases as \L\ increases. 

This remark also leads to the following picture. Having fixed qq = 0, let /o increase 
(or decrease, the sign is physically irrelevant) from /o = 0, and generate configurations for 
{f{u),g{u)) using (^). Since /o sets the value for f{u) (when L = 0) or for its |L|-th order 
derivative (when L 0) at n = 0, and because of the negative driving force —K^f{u) [l — /^(u)] 
(for |/(it)| < 1) in the LG equation, the initial oscillatory pattern of the Bessel function is then 
pushed further and further outwards, with the position of the innermost zero of f{u) moving 
out while the spacing Ar between successive zeros retains more or less the above constant value 
(due to the condensate rigidity set by .^). Consequently, more and more zeros of f{u) leave the 
sample until configurations with no zeros left within it are reached. Moreover, throughout the 

(n) 

increase of /o from the null value, there is a value /q to be found within each of the successive 
intervals in /o which are associated to the loss of a zero in f{u) from the sample, such that the 
b.c. df{ub)/du = in (IT^) is also satisfied. Correspondingly, the solution for g{u) then sets 
the associated value for 6ext through the second b.c. at u = ui, in ([T^). Hence in conclusion, 
with L fixed and go = 0, this demonstrates the existence of (L, n) annular vortex solutions 
with n = 0, 1, • • • ,n^i^x{L), characterized by having the order parameter vanish on concentric 
cylinders when n 7^ and on the symmetry axis when L 7^ 0. 

In the (ftextj^norm) phasc diagram, the above procedure determines a certain curve starting 
from the origin]^. Along that curve, there thus appears in succession a series of points associated 
to the annular vortex solutions with a value of n decreasing from n^g,x{L) to n = 0. Each of 
these points thus also defines a local extremum of the free energy (^) (which may then take 
either a negative or a positive value), and in fact the numerical analysis Q shows that these 
are local minima in the parameter space {fo,9o) (nevertheless, this is a far cry from having 
established dynamical stability for these solutions within the entire configuration space). 



Finally, the complete set of solutions to (|8|) with the b.c. (12) is generated for whatever 
values of 6ext by letting now go run up and down its parameter space starting from go = 0, 
thereby generating energy curves in the (6ext, "fnorm) phase diagram associated to each of the 
{L,n) annular vortex solutions. This is to be done for each of the (/o"\<7o = 0) solutions 
constructed above, by adjusting appropriately the value of /o for each new value of go. To each 
such solution, there corresponds a specific value for 5ext set by the b.c. in (|l2|). However, this 
correspondence between go and 6ext is not necessarily one-to-one, as is indeed demonstrated by 
the numerical analysis. Thus, there may exist more than one solution of fixed (L, n) for some 
values of the external field, a fact responsible for hysteresis and switching phenomena between 
states of different L values in a varying field, associated to the crossing points of the negative 
energy curves of {L,n) vortices in the (6ext; "^norm) phase diagram. 

The Maxwell equation in (^ shows that the magnetic field b(u) is stationnary at the zeros 
of f{u) (there is no Meissner effect at those points) and that it is monotonous in between zeros 
of g{u). Moreover, since g{u) ~ uJ {u) / f"^ {u) is finite at the nodes of /(ti), the supercurrent 
J{u) must also vanish at these locations while it remains of constant sign for as long as g{u) 
does not cross the null value. In particular, when L and 6ext are of opposite sign, g{u) does 
not vanish within the sample, so that all the successive supercurrent vortices retain the same 
circular orientation. 



series expansion analysis of the solutions shows that in the case L = 0, this curve runs exactly along the 
negative £norm axis, in which case go — Q always implies fecxt = 0. 
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The same considerations apply in the case of the annular geometry, the two free b.c. to be 
adjusted in a likewise fashion being then the values of f{u) and g{u) at u = Ua- What replaces 
then the result (^) is a linear combination of the Bessel functions of the first and second kind, 
J|i|(Kti) and N^h{ku), while otherwise all the previous considerations remain applicable, with 
the obvious provision that the total number n of zeros in f{u) is reduced as compared to the 
disk case because of the inner hole of radius • 

Hence, the above analytical considerations, though not providing explicit solutions, demon- 
strate beyond any doubt the existence of the annular vortex solutions with the characteristics 
described in the Introduction. However, the important issue of their dynamical stability {i.e. the 
existence of negative modes of the quadratic operators in (^)) remains open. A remark possibly 
of interest in this respect, to be detailed elsewhere[^, is the following. It is well known[^ that 
for the critical value k = l/\/2, the ordinary giant vortex solutions {L,n = 0) in the infinite 
dimensional plane saturate the Bogomol'nyi bound |P for the free energy by satisfying first-order 
differential equations from which then follows. However, given rather the b.c. (|l2|) and (|l3|) 



associated to a sample of finite radial extent, it may be shown that these first-order differential 
equations are incompatible with these b.c, so that the solutions {L,n ^ 0) cannot be consi- 
dered as being Bogomol'nyi type topological configurations. One may argue that the solutions 
(L,n = 0) are only a slight deformation of the Bogomol'nyi configurations implied by the b.c, 
but this remark does not apply to the {L,n ^ 0) solutions which are certainly not degenerate 
in energy and thus cannot all saturate the same Bogomol'nyi bound. Indeed, as soon as f{u) 
crosses the null value at some radius, from there on at larger radii the function f{u) keeps 
vanishing at almost regular intervals in an oscillatory pattern, so that the free energy of such 
configurations keeps increasing with increasing values of n since each successive supercurrent 
vortex contributes some additional amount of kinetic energy. This is the physical reason why 
annular vortex solutions were never found through formal mathematical analyses of the space 
of solutions to the LG equation in the infinite plane|||]. 

4. Numerical solutions. The above considerations are of course explicitely confirmed through 
the numerical resolution of the equations @ with the b.c. (|T^ and (|l3|). By lack of space, only 
a small — though representative — sample of examples is presented |p in the disk case (those for 
the annulus are very much similar). The method followed is the one described above using the 
free parameters (/o,5o) ^-nd a 4**^ order Runge-Kutta integration of the differential equations. 

Fig.l displays the (^ext, "fnorm) phase diagram for a choice of radius U), = r^/X = 5 and LG 
parameter k = 1. Even though only the curves for L = 0, ±1, ±2, ±3 are shown, the (L, n = 0) 
giant vortex configurations follow of course the usual pattern. Solutions with n = 1 exist only 
for L = 0. Indeed, for L ^ 0, because of the central volume taken up by the magnetic vortex 
at the center (see Fig. 2), the position of the first zero of f{u) is so much pushed outwards that 
it leaves the sample for this choice of radius. For larger values of KUf, = Vb/^, n = 1 solutions 
also exist for L 7^ 0, in which case a series of energy curves similar to those shown in Fig.l 
for the n = solutions also appear for the n = 1 solutions (and so on for still larger values 
of n, provided KUb is large enough; the choice KUb = 5 is made in order that Fig.l still looks 
simple enough to the eye) . Clearly, the crossings of all such curves render the dynamics of such 
superconducting devices of large enough radii still far richer and more complicated than that 
provided already by the usual n = configurations]^. 

Another feature of Fig.l worthwhile to emphasize is the caustic- like structure in the energy 
curves around their crossing points with the zero energy axis. The appearance of these cusps 
stems from the method used to generate configurations. Solutions move along these curves as 
the value of go varies, so that the relation between gQ and 6ext is not necessarily one-to-one for 
some intervals in these quantities. Moreover, solutions then also move into the positive energy 
domain, until they meet a return point back towards the zero energy axis which they reach 
sometimes after having moved back into the negative energy domain as well. Quite clearly, this 



6 



specific feature displayed by our approach[^ should prove to be important to an understanding 
of the hysteresis phenomena which exist around the critical fields for the normal-superconducting 
phase transitions to the {L,n) solutions (even though not apparent on Fig.l, n ^ solutions 
also possess such caustic behaviour), whose physical properties are determined by the dynamic 
and thermodynamic decay rates between the normal and these different superconducting con- 
figurations. These rates are function of the applied field i?ext and the temperature T (note that 
A(r) and ^(T) also depend on the latter parameter), as well as the value and especially the sign 
of the rate of change in the applied field. 

Fig. 2 displays all the configurations which exist with L = 0, 1 for a disk of radius Ub = 11 
and LG parameter k = 1, submitted to an external field of 6ext = 0.05 (for A = 50 nm, this 
corresponds to Sext ~ 66 Gauss). The values of f{u), b{u) and p{u)g{u) /u = q\^J{u)/fL ~ 
— J{u) are shown for the (L = 0; n = 0, 1, 2, 3) and (L = 1; n = 0, 1, 2) solutions. In the latter 
case, the solution with n = 3 has disappeared because of the space taken up by the central 
magnetic vortex, which thus pushes outwards the whole train of zeros of the (L = 0, n = 3) 
solution. Note also that for this disk with u;, = 11 and k = 1, the values of each of the solutions 
n = 0, 1,2 coincide almost exactly beyond u > OAuh when comparing the L = and L = 1 
configurations, since the effects of the L = 1 magnetic vortex remain confined within a radius of 
a few penetration lengths A. Further properties of these configurations have been discussed in 
general terms already, in particular the fact that the magnetic field b{u) is stationnary at nodes 
of the order parameter ^(r, (j)), where the supercurrent J{u) then also vanishes, thereby enabling 
further penetration of the external field within the superconducting sample and a partial anti- 
screening of the Meissner effect for n ^ annular vortices as compared to the usual n = 
configurations. These few examples thus confirm exactly the results advocated on basis of the 
general analytical considerations developed previously. 

5. Conclusions. This letter has demonstrated the existence of annular vortex solutions to 
the Landau-Ginzburg equations in mesoscopic superconductors. These configurations include 
and generalise the well known Abrikosov and giant magnetic vortex ones, are characterized 
by properties described above, and exist only thanks to the finite spatial extent of mesoscopic 
devices. Important issues remain open obviously, most of them raised by the eventual experi- 
mental observation of direct or indirect physical effects related to the existence in actual samples 
of such annular vortices, either as single entities or organised into regular or irregular ensembles. 
Dynamic and thermodynamic stability is obviously essential for their direct observation, and 
remains to be investigated. Irrespective of their stability, these configurations should have indi- 
rect manifestations since they enrich the phase diagram of mesoscopic superconducting disks and 
loops in an external magnetic field in ways still to be explored. Besides hysteresis phenomena in 
such systems to which our approach should be applicable, annular vortices also provide for a new 
mechanism for the switching between states of different fluxoid quantum number, which — since 
these solutions manifestly share the cylindrical symmetry of the sample — is in direct competition 
with the specific mechanism unravelled in Ref.[10U. Static and dynamic properties of mesoscopic 
superconducting devices should thus be affected by the existence of such annular vortices, in 
ways to be assessed explicitely and whose conclusions could potentially be of importance when 
having specific practical applications of such devices in mind. We plan to report on progress on 
some of these issues in later work. 
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Figure Captions 



Figure 1: The normalised free energy ^norm of ( ll^ ) a. function of the normalised applied field 
6ext for a disk with = 5 and k = 1. Displayed are the n = giant vortex configurations with 
L = —3, —2, —1, 0, 1, 2, 3 from left to right in that order (thin parabolic curves), while the n = 1 
solution exists for L = only (thick line). 



Figure 2: The (L = 0; n = 0,1,2,3) and {L = l;n = 0,1,2) configurations for a disk with 
Uh = 11, K = 1 and 6cxt = 0.05. Displayed in that order from left to right are the values for 
/(n), b{u) and f'^{u)g{u) /u = q\^J{u)/h ~ —J{u) as functions of the variable x = u/uf, = r/ri,, 
< X < 1. The top (resp. bottom) panels correspond to L = (resp. L = 1). The /(n) 
values for the (L = 0, n = 0) solution coincide almost exactly with f{u) = 1 and cannot be 
distinguished from unity on these graphs (one has f{uf,) = 0.9995072 and f{uf,) = 0.9995076 for 
(L = 0,n = 0) and (L = l,n = 0), respectively). 
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